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ABSTRACT 

Detection of Jupiter mass companions to nearby solar type stars with precise radial velocity measurements 
is now routine, and Doppler surveys are moving towards lower velocity amplitudes. The detection of several 
Neptune-mass planets with orbital periods less than a week has been reported. The drive toward the search 
for close-in Earth-mass planets is on the agenda. Successful detection or meaningful upper limits will place 
important constraints on the process of planet formation. In this paper, we quantify the statistics of detection 
of low-mass planets in-close orbits, showing how the detection threshold depends on the number and timing 
of the observations. In particular, we consider the case of a low-mass planet close to but not on the 2:1 
mean motion resonance with a hot jupiter This scenario is a likely product of the core-accretion hypothesis 
for planet formation coupled with migration of jupiters in the protoplanetary disk. It is also advantageous 
for detection because the orbital period is well-constrained. We show that the minimum detectable mass is 
« 4 Me (A^/20)-'/2(cr/m ?r^){P /dyl^iM^/Mofl^ for > 20, where is the number of observations, P the 
orbital period, a the quadrature sum of Doppler velocity measurement errors and stellar jitter, and Mi, the stellar 
mass. Detection of few Earth mass rocky cores will require ~ 1 m s"' velocity precision, and most important, 
a better understanding of stellar radial velocity "jitter". 
Subject headings: planetary systems 



1. INTRODUCTION 

Over 100 Jupiter-mass planets have been discovered around 
nearby stars with precise radial velocity measurements (see 
Marcy et al. 2003 for a review). As time goes by, the range of 
periods and amplitudes accessible to these surveys increases, 
moving to longer periods and lower amplitudes (e.g. Carter 
et al. 2003; Fischer et al. 2003; Jones et al. 2003). Recently, 
the detection of Neptune-mass planets with periods less than 
ten days has been reported (Butler et al. 2004; McArthur et 
al. 2004; Santos et al. 2004). In this paper, we quantify the 
detectability of low mass planets in close orbits, and discuss 
the prospects for detecting earth mass objects. We first dis- 
cuss possible origins of such planets and the importance of 
detecting them. 

1.1. The origin and importance of low mass rocky cores in 
close orbits 

In the conventional planet formation scenario (Safronov 
1969; Pollack et al. 1996), heavy elements coagulate into 
terrestrial-planet-like cores prior to the formation of Jupiter- 
mass planets by subsequent accretion of gas. The critical mass 
which segregates these two populations is Mc ~ a few M®, 
determined by the requirement of sufficiently high cooling ef- 
ficiency in the gas envelope (Stevenson 1982; Bodenheimer & 
Pollack 1986; Ikoma et al. 2000). In a minimum mass nebula, 
both type of planets are thought to form at distances from their 
host stars comparable to those of planets in the solar system. 

As a consequence of their tidal interaction, gas giant planets 
may undergo orbital migration (Goldreich & Tremaine 1980; 
Lin & Papaloizou 1986) which is finally halted either by their 
interaction with their host stars or through the local or global 
depletion of disks (Lin et al. 1996; Trilling et al. 2002; Ar- 
mitage et al. 2002; Ida & Lin 2004). In this process, many 
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gas giants are expected to be disrupted, transferring their mass 
onto the host star (Trilling et al. 2002), perhaps driven by tidal 
inflation (Gu, Lin, & Bodenheimer 2003). Ida & Lin (2004) 
estimate that 90-95% of planets that migrate to a < 0.05 AU 
must perish. Recent models of evaporation of material from 
hot jupiters (Lammer et al. 2003; Baraffe et al. 2004; Lecav- 
alier des Etangs et al. 2004) find that gas giants may undergo 
significant mass loss, and perhaps be completely evaporated, 
during their lifetimes, although the efficiency and rates of 
evaporation remain uncertain (Yelle 2004). This inference 
raises the possibility of a class of remnant "hot Neptunes" or 
rocky cores in close orbits^. 

An alternative way to produce low-mass planets in close or- 
bits is by resonant capture during migration of a hot jupiter 
Along their migration paths, gas giant planets capture other 
planets onto their mean motion resonances as in the case of 
GJ 876 (Lee & Peale 2002; Nelson & Papaloizou 2002; Kley 
et al. 2004). Mandell & Sigurdsson (2003) also consider the 
survivability of earth mass objects in the wake of Jupiter mi- 
gration. A natural implication of the core-accretion followed 
by migration scenario for the origin of hot jupiters is that they 
may harbor earth-mass planets ("hot earths") orbiting close to 
their mean motion resonances (Aarseth & Lin 2004, in prepa- 
ration). In the gravitational instability scenario, however, both 
ice giants and terrestrial planets are assumed to emerge long 
after the formation of gas giant planets (Boss, Wetherill, & 
Haghighipour 2002). 

Due to the effect of tidal circularization, the semimajor axis 
of the hot earth is expected to be slightly interior to the mean 
motion resonance, and its survival is ensured by the relativis- 
tic precession induced by the gravitational potential of the 
host star (Novak et al. 2003, MardHng & Lin 2004). After 
halting its inward migration, a hot jupiter may move outwards 
due to either its tidally induced mass loss (Trilling et al. 2002; 

^ In this paper, we use the adjective "hot" to describe a closely-orbiting 
planet (orbital period < 10 days). 
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Gu et al. 2003), or tidal interaction with a rapidly spinning 
host star (Stassun et al. 1999; Dobbs-Dixon et al. 2004). In 
both cases, the periods of the close-in earth and hot jupiter 
may not be close to commensurability. 

The presence of earth-mass planets close to the external 
mean motion resonance of a hot jupiter is also possible. In the 
wake of a migrating gas giant, the prototerrestrial planets in 
the disk external to its orbit may also undergo type I inward 
migration (Ward 1986) as a consequence of their own tidal 
interaction with the disk gas (Goldreich & Tremaine 1980). 
By inducing the formation of a gap near its orbit, a migrat- 
ing jupiter provides a tidal barrier which prevents any migrat- 
ing terrestrial planets from passing it. But this possibility is 
uncertain since type I migration is strongly affected by the 
poorly known turbulence near the coorbital region (Koller et 
al. 2003; Nelson & Papaloizou 2004). If this process is effi- 
cient, it could lead to the formation of an isolated hot earth 
(Ward 1997; Bodenheimer et al. 2000). Even if type I mi- 
gration is inefficient, the presence of a hoi jupiter can induce 
planetesimals to accumulate and form a hot earth just beyond 
the outer edge of the gap (Bryden et al. 2000). In the pres- 
ence of additional long-period gas giant planets, a sweeping 
resonance may also lead to the inward migration and accu- 
mulation of terrestrial planets just exterior to the orbits of hot 
Jupiters (Lin, Nagasawa, & Thonraies 2004). 

Based on the core-accretion scenario for planet formation, 
it is already tempting to infer the existence of terrestrial plan- 
ets in extrasolar systems from the detection of jupiter-mass 
gas giants (Ida & Lin 2004). However, the actual detection 
of a terrestrial mass object would be the first confirmation of 
rocky cores in an extrasolar planetary system (see however 
Kuchner 2003). The association of a hot earth with the in- 
terior mean motion resonance of a hot jupiter is particularly 
interesting because it would indicate that cores are formed 
prior to the emergence and migration of gas giant planets, pro- 
viding strong support for the conventional core-accretion sce- 
nario for planet formation, and for the ubiquity of terrestrial 
planets in extrasolar planetary systems. The presence of a hot 
earth exterior to the orbit of a hot jupiter would support the 
notion of a tidal barrier which has the potential to enhance the 
formation probability of multiple planet systems. Finally, the 
discovery of any isolated hot earth would support the concept 
of type 1 migration. 

1.2. Outline of this paper 

In this paper, we discuss the detectability of such low mass 
companions in radial velocity surveys. In particular, we con- 
centrate on the scenario in which a terrestrial planet was cap- 
tured into the mean motion resonance of a migrating jupiter- 
mass planet, resulting in the formation of a "hot earth" which 
is near the interior mean motion resonances of a "hot jupiter". 
With the limitations of radial velocity surveys, the detection 
of a hot earth external to the orbit of a hot jupiter is much 
more difficult. We also calculate the detection thresholds for 
isolated low mass planets. 

Being offset from the mean motion resonance of any nearby 
hot Jupiters, the expected contribution to the radial velocity 
amplitude K from a closely orbiting several Earth-mass object 
is 
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where P is the orbital period, Mp is the mass of the planet, 
and is the mass of the star. This is comparable to the 
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Fig. 1. — Simulated velocity residuals (top panel) after subtraction of a 
0.5 Mj mass planet with P = 3.1 days, K = 70 m/s. We simulate one observa- 
tion per night for 30 nights, with measurement error of 1ms"'. Periodogram 
(lower panel) evaluated for periods from 1 to 4 days before (dotted hne) and 
after (solid line) subtraction. The strong feature present at 1.5 days before 
subtraction is due to aliasing of the 3.1 day signal. Subtracting reveals the 
presence of a 3 companion at P = 1.4 days. The dashed line shows the 
detection threshold. 



expected velocity variability from measurement errors, and 
intrinsic stellar "jitter" (Saar, Butler, & Marcy 1998). The 
orbital frequency is well-predicted for a hot earth close to a 
mean motion resonance with a hot jupiter, and so we expect 
the detectability to be somewhat better than a blind frequency 
search. Here we calculate the detectabihty and discuss the 
number of observations, and observing strategy needed for 
detection. 

2. DETECTABILITY OF CLOSE COMPANIONS 

2.1. Search Technique 

We simulate a set of observational data for a system of plan- 
ets which are on circular orbits. We search for periodicities 
using the Lomb-Scargle periodogram (Lomb 1976; Scargle 
1982). For a set of A'^ observation times {tj}, velocities {v^}, 
and measurement errors {aj}, and a trial orbital frequency 
u} = 2Tr/P,we fit the function 



fj = A cos uitj + B sin cotj + C 



(2) 



to the velocity data by minimising = ( 1 / jy) ^j(vj-fj)^/(Tj. 
The number of degrees of freedom is iy = N-3 since there are 
3 parameters in the model. We measure the goodness of fit as 
a function of trial frequency co using the periodogram power 
z, defined as 

AxV2 



(3) 



where Ax^ = (A^- 1)Xa'-i ~ '^X?" ^^'^ Xn-\ is the reduced of 
a fit of a constant to the data, (A^- 1)Xa'-i = X)/^/ ~ M^'/'^j- 
Here, we extend the original Lomb-Scargle periodogram by 
allowing the mean to float at each frequency (Walker et 
al. 1995; Nelson & Angel 1998; Gumming, Marcy, & Butler 
1999), rather than subtracting the mean of the data (v) prior to 
the fit. The periodogram power measures the improvement of 
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when the sinusoid is included in the fit, in a similar way to 
a classical F-test (see Gumming 2004 for a recent discussion). 

A typical Doppler measurement error is 3-5 m/s, although 
the precision of Doppler surveys continues to improve to- 
wards 1 m/s (Mayor et al. 2003; Butler et al. 2004). How- 
ever, at these high precisions, stellar "jitter" becomes the lim- 
iting factor Jitter is observed at the few m/s level, depend- 
ing on stellar properties such as age, rotation, and level of 
magnetic activity (Saar et al. 1998; Santos et al. 2000). For 
simplicity, we add Gaussian noise with amplitude a = 1 m/s 
throughout this paper, which represents both Doppler errors 
and stellar jitter. The mass which can be detected scales pro- 
portional to a, so that our results can be rescaled to the appro- 
priate value of a. For example, a 5 planet with ct = 5 m/s 
has the same detectability as a 1 planet with a = 1 m/s. 
We have checked that this scaling applies for the two planet 
case, when the mass of the hot jupiter is much greater than the 
hot earth. 

Our search procedure is to find the best-fit orbit for the 
hot jupiter, subtract it from the data, and search the residu- 
als for periodicities. This is a much simpler approach than 
simultaneous fitting of both sinusoids, and is valid as long as 
the parameters from the two fits are uncorrelated. We have 
checked that this is the case, and find that only for sampling 
rates much greater than 1 day"' is there some improvement 
by simultaneous fitting. An alternative possibility is that the 
orbital solution for the hot jupiter is already known from pre- 
vious measurements, in which case this information can be 
used to subtract the hot jupiter's signal. Therefore, we also 
consider the cases where the full orbital solution or only the 
orbital period is known in advance. 

If the hot jupiter orbital period is known in advance, we 
find the best fit amplitude and phase at that period. If no in- 
formation is available, we search for the hot jupiter by eval- 
uating zipj) for periods between 2 and 4 days on a grid with 
frequency spacing l/4r, where T is the duration of the ob- 
servations. To adequately subtract the hot jupiter's signal, 
it is important to accurately determine its orbital frequency. 
Therefore the frequency with the largest periodogram power 
is used as the starting value for a more accurate search for the 
maximum of z{lo). Once the hot jupiter's orbit has been sub- 
tracted, we calculate z((jj) for the residuals, and find the maxi- 
mum periodogram power z^ax near the 2: 1 resonance. We do 
not include a constant term in the fit to the residuals, giving 
v = N-2. Figure[2shows an example of the velocity residuals 
and the periodograms before and after subtraction. The peak 
at w 1.5 days in the initial periodogram is due to aliasing of 
the hot jupiter's frequency. After subtraction, the signal due 
to a 3 planet at 1 .4 days shows clearly as a peak in z- 

The importance of accurately determining the hot jupiter's 
orbital period is illustrated by Figure|2 which shows two ex- 
amples of the velocity residuals after subtraction of the hot 
jupiter signal. In the first case (upper panels), we take the hot 
jupiter period returned by the periodogram (evaluated with a 
frequency spacing 1 /AT); in the second case (lower panels), 
we find a more accurate period estimate by searching for the 
peak in z{lo). In the first case, the velocity residuals are dom- 
inated by the inadequately-subtracted hot jupiter signal rather 
than the hot earth. How accurately must the hot jupiter period 
be known? If the orbital frequency is known to an accuracy 
5uo, the amplitude of the residual part of the hot jupiter's sig- 
nal is Av A\5loT, where A\ is the hot jupiter amplitude. 
The accuracy to which the hot jupiter frequency may be de- 



termined is 5uj « (27r/r)(CT/%/A^Ai) (e.g. Bretthorst 1988), so 
that 



independent of T and A i . For large A^, Av <C ct, so that the 
residual part of the hot jupiter signal has no effect on de- 
tectability. However, for < (27r)^, the residuals from the 
subtraction contribute a significant additional source of ve- 
locity variability. Even when the hot jupiter period is known 
in advance, there is an uncertainty in the fitted amplitude A\ 
of w g/\/N, and phase of w (2/Ny^^(a'/K), giving an addi- 
tional velocity scatter for low A^. Therefore for small we 
expect detection of a hot earth to be more difficult when the 
full orbital solution for the hot jupiter is not specified in ad- 
vance. 

2.2. Calculation of Detection Threshold and Detection 
Probabilities 

The significance of the maximum observed periodogram 
power Zmax depends on how often an equally good or better fit 
would occur purely due to a noise fluctuation. We determine 
this with Monte Garlo simulations. We generate data sets with 
a hot jupiter plus noise, i.e. without a hot earth, and search for 
a second companion. The 99% detection threshold Zd is deter- 
mined as the value of z which is exceeded in only 1 % of trials, 
or alternatively for which there is a 1 % false alarm probabil- 
ity F. The detection threshold is indicated in Figure^by the 
dashed line. 

We adopt a similar Monte Garlo approach to determine the 
detection probability. We generate a large number of data sets 
for each choice of Mp and A^, and calculate the fraction of tri- 
als for which the hot earth is detected. A random distribution 
of inclinations is included. Figure |3l shows the detection effi- 
ciency as a function of number of observations A^ for Mp = 1, 
1.5, 2, and 3 M®. We simulate 1 observation per 8 hour night, 
and take observations for successive nights. In these simula- 
tions, we search for the hot earth orbital period between 1 and 
1.8 days. Figurel^shows the effect of changing this frequency 
range. If the frequency of the hot earth is specified in ad- 
vance, the detectability is increased since the number of "in- 
dependent frequencies" in the search is less, giving a smaller 
chance of a false alarm due to a noise fluctuation. This is the 
well-known "bandwidth penalty" (e.g. Vaughan et al. 1994). 

In Figure|5l we show the number of observations needed to 
detect a given mass 50% of the time. The crosses are for data 
sets containing only a hot Earth, or equivalently, for the case 
when the hot jupiter orbit is completely known in advance. 
The triangles are for the case with a hot jupiter whose or- 
bital period is known in advance, and the circles for the case 
with a hot jupiter, but with a full search for both orbital pe- 
riods. For A^ « 10-20, planets with M > 4(ct/1 m s"') 
can be detected. However, detection of a planet with M ^ 
1 (ct/1 m s"') requires A^ w 200. For A^ < 20, detection of 
the hot earth after fitting for the hot jupiter's signal is harder 
than detection of the hot earth alone or with the hot jupiter or- 
bit fully specified in advance (compare the crosses and circles 
in Fig|5}- This is likely due to the additional velocity scatter at 
low A^ from inadequate subtraction of the hot jupiter's signal 
(Av > CT in eq. El). 

2.3. Analytic estimates 

In this section, we derive an analytic estimate for the detec- 
tion threshold, following the approach of Gumming (2004), 
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Fig. 2. — The velocity residuals and periodogram after subtraction of the hot jupiter for two different estimates of the hot Jupiter's orbital frequency. In the 
upper panel, we use the frequency estimate from the periodogram evaluated with frequency spacing 1/47. In this case, the poor subtraction leads to a large scatter 
in the residuals. In the lower panel, we determine an accurate estimate of the frequency by finding the peak of z{uj). The velocity residuals are then dominated by 
the 5 M© companion. 
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Fig. 3. — Detection efficiency for hot earths with masses of I, 1.5, 2, and 
3 M0 and period 1.4 days, with cr = 1 m s"'. We assume 1 randomly-timed 
observation on successive 8 hour nights. We use 10,000 trials to evaluate the 
detection probability. The hot Jupiter has a mass of 1.0 Jupiter mass, and a 
period of 3.0 days. 
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Fig. 4. — The effect of the range in orbital period searched on the detection 
probability. We take a = 1ms"' and Mp = 3 Mg. Triangles represent orbital 
periods ranging from 1-4 days, squares 1.3 - 1.5 days, and circles 1.39 - 
1.41 days. The dotted line is for 1-1.8 days, which is used elsewhere in this 
paper. As a larger frequency range is searched, iV,- increases giving a larger 
detection threshold and lower detection probabUity. 



who discusses the detectability of single planets with radial 
velocities. First, we determine Zd semi-analytically for Gaus- 
sian noise. The cumulative distribution of for a single 
frequency is (e.g. Cumming et al. 1999) 

Prob(z>zo)= + ^ , (5) 

or for large N, 

Prob(z > zo) ~ exp(-zo), . (6) 



For a given false alarm probabihty F, the detection threshold 
is given by 

F=l-[l-Prob(z>Zd)fS (7) 

where A^, is the number of "independent frequencies" 
searched. For F <^1, 

F^NiPwhiz>Zd). (8) 

For unevenly-sampled data, A'; must be determined by Monte 

Carlo simulations. However, since the spacing of peri- 
odogram peaks is 1/r, a rough estimate is A^, « TAf, where 
A/ is the frequency range searched (Cumming 2004). 
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Fig. 5. — The number of observations needed to detect a planet 50% of the 
time as a function of the mass Mp scaled by the en'or cr. Observations are 
taken at a random time during each consecutive night. The crosses are for 
data sets containing only a hot Earth, the triangles for the case where the hot 
jupiter period is known in advance, and the circles for a full search for both 
orbital periods. The curves are the analytic results. The dotted curve occurs 
if the number of independent frequencies, Nj, is Nj = TAf, where T is the 
duration and A/ is the frequency range used to search for the hot earth. For 
the dashed curve, W, = 10. 



We next make an analytic estimate of the detection thresh- 
old. In the presence of a signal with amplitude K, the average 
power is z, = iyK^/4a^ (Groth 1975; Scai-gle 1982; Home & 
Baliunas 1986, we have accounted for the different normal- 
ization). For a given A^,- and F, we find the detection threshold 
Zd by inverting equation ([7) using the analytic distribution (|5}- 
Then, setting z, = Zd gives the signal to noise ratio needed to 
detect the signal 



K 
V2a 



2/(iV-3) 



21n(A^,/F)"l'/^ 



1/2 



N 



(9) 
(10) 



(Gumming et al. 2003; Gumming 2004) where the approxi- 
mation is for large A^. As expected, the signal to noise ratio 
goes down as 1 / \/N for 1 . 

Gombining equation ([0 with equation (|9} or (I10> gives the 
velocity amplitude and mass which can be detected 50% of 
the time. Equation (I10> (A^ ^ 1) gives 
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fHNi/F] 
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lOMff 
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Me 



(11) 



2/3 



(12) 

where we have also used the fact that the mean value of sin/ 
is 7r/4, and we take A^,- = 100 and F = 0.01. 

Figure [S] compares the analytic results for arbitrary A^ 
(eq. |9l) with our numerical simulations. The dashed and dot- 
ted curves show the analytic models for two different esti- 
mates of A^,. The agreement is excellent for a single planet. 



For the two planet case, equation ( I12t applies for A^ > 20, but 
underestimates M50 for A^ < 20. For such low A^, we expect 
Av > a (eq. 0), making detection of the hot earth more dif- 
ficult than in the single planet case. 

2.4. Observing Strategy 

We now address to what extent the spacing of the observa- 
tion times affects the detectability of the hot earth. First, we 
consider observations made on successive nights, but with dif- 
ferent numbers of observations per night. Figure|6lshows the 
detection probability as a function of number of nights, and as 
a function of number of observations. For a fixed number of 
observations, increasing the sampling rate reduces the dura- 
tion of the data, leading to less frequency resolution and there- 
fore smaller A^, . We therefore expect better detectability for a 
faster sampling rate at a fixed number of observations. This 
is initially the case, as can be seen by comparing the curves 
for 1 per night and 5 per night in Figure |6l However, for 
very rapid sampling, a large number of observations must be 
made before a complete hot jupiter orbit is sampled, therefore 
allowing the orbital parameters to be determined accurately 
enough for a good subtraction. This is likely the reason for 
the decrease in detectability on going to more than 5 observa- 
tions per night in Figure|6l Although observing more rapidly 
always leads to detection in a fewer number of nights (left 
panel of Fig. |6j, the number of observations required even- 
tiually becomes very large (right panel of Fig.|6jl. Therefore 
there is an optimum observing rate (roughly a few per night 
for this case, but it depends slightly on the hot earth mass). 

Although observing on successive nights would be possi- 
ble with a dedicated telescope, current radial velocity surveys 
are limited by telescope scheduling. Therefore, we have also 
simulated more realistic observation times for current surveys, 
by observing for three nights each month, with the seperation 
between observing runs randomly chosen between 20 and 40 
days. We find that if the hot jupiter orbital parameters are 
not known in advance, the detectability can be affected by the 
observing scheme, because of aliasing of the hot jupiter fre- 
quency leading to subtraction of an incorrect orbit. However, 
in the practical case that the hot jupiter period is known from 
previous observations, the hot jupiter orbit can be adequately 
subtracted, and the effect of the observing strategy is small. 

3. SUMMARY AND DISCUSSION 

We have calculated the detectability of low mass planets in 
radial velocity surveys, and in particular a "hot earth" com- 
panion to a hot jupiter. Detection of such a companion would 
give important clues to the ordering of the planet formation 
process, and the ubiquity of terrestrial mass cores around so- 
lar type stars. The velocity amplitude and mass required for 
a 50% detection rate are given by equations il 1> and ( I12t . 
and shown in Figure |5] For A^ « 20, masses greater than 
4 Me (cr/m s-^)(P/dy^\M^/MQff^ have 50% detectability 
or better 

Our results apply to a low mass planet that is near the mean 
motion resonance with a hot jupiter, or that is isolated. In 
Figure we show the detection limits of equation J12t com- 
pared to the present distribution of known exoplanets"*. The 
left panel is for = 1 Mq, the right panel for Mj, = 0.2 Mq. 
Lower mass stars have a greater velocity amplitude for a given 
planet mass, so that M dwarfs are the most promising to 

Taken from http://www.exoplanets.org/. 
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Fig. 6. — Detection probability for a 1 planet with different numbers of observations made on consecutive nights. We show the detection probability 
against number of nights (left panel) and number of observations (right panel). A more rapid sampling rate leads to a detection in fewer nights, but at the cost of 
taking many observations. 



TABLE 1 

Summary of radial velocity data for hot iupiters and estimated upper limits for hot earths 



Star name 


References 




rms'' 


N 


T 


W 




Kg/ 








(d) 


(ms-') 




(d) 




(Mq) 


(m s-') 


(M(b) 


HD 73256 


1 


2.55 


15 


40 


80 


4 


1.05 


22 


58 


HD 83443 


2,3 


2.99 


3.8 


36 


1150 


38 


0.79 


7.3 


15 


HD 46375 


4 


3.02 


2.6 


24 


516 


21 


1.0 


6.5 


15 


HD 179949 


5 


3.09 


10 


23 


735 


30 


1.24 


27 


74 


HD 187123 


6 


3.10 


7.6 


20 


250 


10 


1.0 


20 


48 


T Boo 


7 


3.31 


16.5 


58 


4400 


166 


1.3 


26 


74 


BD- 103 166 


8 


3.49 


8.1 


17 
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"Orbital period of the hot jupiter. 

''rms of the residuals to the orbital solution. 

'^Estimated number of independent frequencies for a search from 2-2.4 times the hot jupiter orbital frequency. A'; ~ O.IT/Pj. 

"•Estimated 99% upper limit, using equation j5) with Kgg « 1 JK^q . 

■^99% upper limit on the mass of a hot earth near the 2: 1 resonance with the hot jupiter. 

'References.— (1) Udry et al. 2003; (2) Mayor et al. 2002; (3) Butler et al. 2002; (4) Marcy et al. 2000; (5) Tinney et al. 2001; (6) Butler et al. 1998; (7) Fischer 
et al. 2001; (8) Buder et al. 2000; (9) Udry et al. 2000; (10) Tinney et al. 2003 



searcli for terrestrial mass bodies. Close orbits are particu- 
larly interesting for this case, since they lie within the hab- 
itable zone for M dwarfs (orbital periods of a few to tens of 
days; Kasting, Whitmire, & Reynolds 1993; Joshi, Haberle, 
& Reynolds 1997). 

Recently, three Neptune-mass planets have been discov- 
ered in close orbits around GJ 436 (Butler et al. 2004; 
Msin/ = 21 Me, P = 2.644d, = 42, rms=5.3 m s"'), p 
Cnc (McArthur et al. 2004; Msin; = 14.2 M®, P = 2.808d, 
N = 119, rms=5.4 ms"'), and p Ara (Santos et al. 2004; 
Msinj = 14 Me, P = 9.5d, = 24, rms=0.9 m s"'). These 
can be seen as the lowest mass planets in Figure0 Only one 
of these stars, GJ 436, is an M dwarf. All of these detections 
lie above the detection limits in Figure^ 

The lowest mass planet detected prior to these recent 
announcements was HD 49674b (Butler et al. 2002) with 
M sin / = 0. 12 My (or w 40 Me), = 1 3 m S-' , and P = 5 days. 



The velocity curve presented by Butler et al. (2002) has 
= 24, and the residual rms is w 5 m s"'. With these values 
of and cr, the thresholds we derive here imply that planets 
with half the mass of HD 49674b (« 10-20 Me depending 
on orbital period) should be detectable, consistent with the 
recent detections. Endl et al. (2003) report 22 velocity mea- 
surements of the M dwarf Proxima Cen with Doppler errors 
of 2.5 m s"'. They calculate an upper limit of 4-6 Me for 
planets in the habitable zone (estimated to be 4-14 days), in 
good agreement with Figure 

A summary of published data for hot jupiters with orbital 
periods Pj < 4- days is shown in Table 1 . To estimate an up- 
per limit on the velocity amplitude of a hot earth, we first use 
the duration and number of observations to estimate A^, for 
a search between orbital periods P//2.4 and Pj/2. We then 
write the 99% upper limit as Ki)<) = 1 .IK^o, where Kso is given 
by equation (|9}, but with the rms of the residuals to the best- 
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Fig. 7. — Summary of detection thresholds in the mass-period plane. The dashed lines show 50% detection thresholds for different N and a, for M* = 1 Mq 
(left panel) and M* = 0.2 Mq (right panel). It is assumed that the duration of the observations is longer than the orbital period. The circles show currently 
detected planets (Msini) from exoplanets.org, and including the 3 recently announced Neptune-mass candidates. The dotted lines show velocity amplitudes of 
K= 1,3, and 10 m/s (sin; =1), and the dashed line shows an approximate detection threshold for SIM (assuming 1/xas sensitivity and a 10 pc distance; see Ford 
& Tremaine 2003). 



fitting hot jupiter orbit substituted for a. Since we adopt the 
full rms of the residuals as the noise level, these are conserva- 
tive upper limits. The values of Kgt) range from 7-30 m s"', 
corresponding to w 15-80 for companions close to 2:1 
resonance with the hot jupiter HD 83443 and HD 46375 have 
upper limits of 15 because of very precise measurements 
giving a low rms (4 and 3 m s"' respectively). 

These detections and upper limits imply that the current 
threshold for detecting planets in close orbits is > 1O-2OM0, 
with planets close to this threshold just starting to be discov- 
ered. With more observations, and improvements in measure- 
ment precision, this threshold can be pushed lower In Figure 
we show the detectable mass for =100 and a = 3 m s~', 
which reaches a few Earth masses for orbits ^ 1 day for a solar 
mass star, or 1 for a 0.2 Mq M dwarf. However, a better 
understanding of stellar jitter will be crucial for convincing 
detections of small amplitude planets. Detection of extrasolar 
planets with radial velocities becomes more and more difficult 
as K approaches a, because alternative explanations for the 
detected signal must be excluded (see the discussion in Gum- 
ming et al. 1999 for example). These include systematic er- 
rors in the measurements, as well as stellar jitter Both may be 
periodic, with timescales associated with seasonal variations, 
or stellar properties such as rotation, convective motions, or 
the appearance and disappearance of magnetic features. Saar 
& Fischer (2000), Paulson et al. (2002), and Saar (2003) dis- 
cuss attempts to correct for stellar jitter using simultaneous 
measurements of activity indicators or spectral line shapes, 
but such studies are just beginning. Very precise Doppler 
measurements with errors at the 1 m s"' level will help to un- 
derstand the stellar jitter and systematic errors. 

Finally, we have assumed that both planets are in circular 
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orbits in this paper. For orbital periods as short as a few days, 
the timescale for tidal circularization is expected to be much 
shorter than the age of the system. However, additional plan- 
ets in the system, at longer orbital periods and so far unde- 
tected, may excite the eccentricity of the hot jupiter Tidal 
damping of an eccentricity induced in this way has been sug- 
gested as the reason for the inflated radius of the transiting 
planet HD 209458b (Bodenheimer, Lin, & Mai'dling 2001). 
The required eccentricity is comparable to current measure- 
ment uncertainties, typically 6e ^ 0.03. Fitting a circular or- 
bit to a Keplerian orbit with eccentricity e gives a residual 
scatter w eAi, where Ai is the amplitude. For a hot jupiter 
with Al w 100 m s"' the scatter in the residuals is compara- 
ble to the signal from a hot earth. However, since the Fourier 
components of a Keplerian orbit are at the orbital period and 
its harmonics, whereas the hot earth is expected to lie out- 
side the 2:1 resonance, it should be possible to distinguish 
between these two possibilities, although simultaneous fitting 
of the hot jupiter and earth orbits may be required. Therefore 
we do not expect a significant reduction in sensitivity over 
the detection thresholds calculated in this paper. This issue 
is closely related to the important question of determining the 
uncertainty in orbital parameters (Ford 2003), particularly the 
eccentricity of hot jupiters, and deserves further study. 



We acknowledge support from NSF through grant AST- 
9987417 and NASA through grant NAG5-13177. AC is sup- 
ported by NASA through Hubble Fellowship grant HF-01 138 
awarded by the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in 
Astronomy, Inc., for NASA, under contract NAS 5-26555. 



Bodenheimer, P & Pollack, J. B. 1986, Icarus, 67, 391 
Bodenheimer, P., Hubickyj, O. Lissauer, J. J. 2000, Icarus, 143, 2 
Bodenheimer, P, Lin, D. N. C, & Mardhng, R. A. 2001, ApJ, 548, 466 
Boss, A. P, Wetherill, G. W., & Haghighipour, N. 2002, Icarus, 156, 291 



8 



Bretthorst, L. G. 1988, "Bayesian Spectrum Analysis and Parameter 
Estimation", Lecture Notes in Statistics vol. 48 (Springer- Verlag) 
Biyden, G., Rzyczka, M., Lin, D. N. C, & Bodenheimer, P 2000, ApJ, 540, 
1091 

Butler, R. P, Marcy, G. W., Vogt, S. S., & Apps, K. 1998, PASP, 110, 1389 
Buder, R. P, Vogt, S. S., Marcy, G. W., Fischer, D. A., Henry, G. W., & Apps, 
K. 2000, ApJ, 545, 504 
Buder, R. P, et al. 2002, ApJ, 578, 565 

Buder, R. P, Bedding, T. R., Kjeldsen, H., McCarthy, C, O'Toole, S. J., 
Tinney, C. G., Marcy, G. W., & Wright, J. T. 2004, ApJ, 600, L75 
Buder, R. P, Vogt, S. S., Marcy, G. W., Fischer, D. A., Wright, J. T., 
Henry, G. W., Laughlin, G., & Lissauer, J. J. 2004, ApJ, submitted 
( astro-ph/0408587 i 

Carter, B. D., Buder, R. P, Tinney, C. G., Jones, H. R. A., Marcy, G. W., 
McCarthy, C, Fischer, D. A., & Penny, A. J. 2003, ApJ, 593, L43 
Gumming, A. 2004, MNRAS, in press i astro-ph/0408470 1 
Gumming, A., Marcy, G. W., & Buder, R. P 1999, ApJ, 526, 890 
Gumming, A., Marcy, G. W., Buder, R. R, & Vogt, S. S. 2003, in ASP 
Conf. Ser. 294, Scientific Frontiers in Research on Extrasolar Planets, 
ed. D. Deming, & S. Seager (San Francisco: ASP), 27 
Dobbs-Dixon, 1., Lin, D.N.C., & Mardling, R.A. 2004, ApJ, in press 
Endl, M., Kurster, M., Rouesnel, F, Els, S., Hatzes, A. P, & Cochran, W. D. 
2003, in ASP Conf. Ser. 294, Scientific Frontiers in Research on Extrasolar 
Planets, ed. D. Deming, & S. Seager (San Francisco: ASP), 75 
Fischer, D. A., Marcy, G. W., Butler, R. P, Vogt, S. S., Frink, S., & Apps, K. 
2001, ApL 551, 1107 

Fischer, D. A., Butler, R. P, Marcy, G. W., Vogt, S. S., & Hemy, G. W. 2003, 
ApJ, 590, 1081 

Ford, E. B. 2003, preprint ( astro-ph/030541 II 
Ford, E. B. & Tremaine, S. 2003, PASP, 115, 1171 
Goldreich, P, & Tremaine, S. 1980, ApJ, 241, 425 
Groth, E. J. 1975, ApJS, 29, 285 

Gu, P, Lin, D. N. C, & Bodenheimer, P H. 2003, ApJ, 588, 509 

Home, J. H. & Baliunas, S. L. 1986, ApJ, 302, 757 

Ida, S. & Lin, D. N. C. 2004, ApJ, 604, 388 

Ikoma, M., Nakazawa, K. & Emori, E. 2000, ApJ, 537, 1013 

Jones, H. R. A., Butler, R. P, Tinney, C. G., Marcy, G. W., Penny, A. J., 

McCarthy, C, & Carter, B. D. 2003, MNRAS, 341, 948 

Joshi, M. M., Haberle, R. M., Reynolds, R. T. 1997, Icarus, 129, 450 

Kasdng, J. F, Whitmire, D. P, & Reynolds, R. T. 1993, Icarus, 101, 108 

Kley, W., Peitz, J. & Bryden, G. 2004, A&A, 414, 735 

Roller, J., Li, H., & Lin, D.N.C. 2003, ApJ, 596, 91 

Kuchner, M. J. 2003, ApJ, 596, L105 

Lammer, H., Selsis, F, Ribas, 1., Guinan, E. F, Bauer, S. J., & Weiss, W. W. 
2003, ApL 598, L121 

Lecavalier des Etangs, A., Vidal-Madjar, A., McConnell, J. C, & Hebrard, 

G. 2004, A&A, 418, LI 

Lee, M.H. & Peale, S. 2002, ApJ, 567 596 

Lin, D. N. C, Bodenheimer, P & Richardson, D. 1996, Nature, 380, 606 
Lin, D. N. C, Nagasawa, M., & Thommes, E. 2004, Nature, submitted 
Lin, D. N. C, & Papaloizou, J. C. B. 1986, ApJ, 309, 846 
Lomb, N. R. 1976, Ap&SS, 39, 447 



Mandell, A. M., & Sigurdsson, S. 2003, ApJ, 599, LI 11 

Marcy, G. W., Buder, R. P, & Vogt, S. S. 2000, ApJ, 536, L43 

Marcy, G. W., Buder, R. P, Fischer, D. A., & Vogt, S. S. 2003, in ASP 

Conf. Ser 294, Scientific Frontiers in Research on Extrasolar Planets, 

ed. D. Deming, & S. Seager (San Francisco: ASP), 1 

Mardling, R. & Lin, D.N.C. 2004, ApJ, submitted 

Mayor, M., Naef, D., Pepe, F, Queloz, D., Santos, N. C, Udry, S., & Burnet, 
M. 2002, in ASP Conf Ser., Planetary Systems in the Universe, ed. A. J. 
Penny, P. Artymowicz, A. M. Lagrange, & S. S. Russell (San Francisco: 
ASP), in press 

Mayor, M., et al. 2003, The Messenger, 1 14, 20 

McAi-thur, B. E. et al. 2004, ApJ, in press ( astro-ph/04085S5} 

Nelson, A. F & Angel, J. R. P 1998, ApJ, 500, 940 

Nelson, R. P, & Papaloizou, J. C. B. 2002, MNRAS, 333, 26 

Nelson, R. P, & Papaloizou, J. C. B. 2004, MNRAS, 350, 849 

Novak, G. Lai, D., & Lin, D.N.C. 2003, in ASP Conf Ser. 294, Scientific 

Frontiers in Research on Extrasolar Planets, ed. D. Deming, & S. Seager (San 

Francisco: ASP), 177 

Novak, G. Lai, D., & Lin, D.N.C. 2004, ApJ, submitted 

Paulson, D. B., Saar, S. H., Cochran, W. D., & Hatzes, A. P 2002, ApJ, 124, 

572 

Pollack, J. B., Hubickyj, O., Bodenheimer, P., Lissauer, J. J., Podolak, M., & 
Greenzweig, Y. 1996, Icarus, 124, 62 

Saar, S. H. 2003, in ASP Conf. Ser. 294, Scientific Frontiers in Research on 
Extrasolar Planets, ed. D. Deming, & S. Seager (San Francisco: ASP), 65 
Saar, S. H., Buder, R. P, & Marcy, G. W. 1998, ApJ, 498, L153 
Saar, S. H. & Donahue, R. A. 1997, ApJ, 485, 319 
Saar, S. H. & Fischer, D. 2000, ApJ, 534, L105 

Safronov, V. 1969, Evolution of thr Protoplanetary Cloud and Formation of 
the Earth and Planets (Moscow: Nauka Press) 

Santos, N. C, Mayor, M., Naef, D., Pepe, F, Queloz, D., Udry, S., & Blecha, 
A. 2000, A&A, 361,265 

Santos, N. C, et al. 2004, A&A, submitted tastro-ph/040847It 
Scargle, J. D. 1982, ApJ, 263, 835 
Schwarzenberg-Czemy, A. 1998, Baltic Astron., 7,43 
Stevenson, D. J. 1982, P&SS, 30, 755 

Stassun, K. G., Mathieu, R. D., Mazeh, T, & Vrba, F J. 1999, AJ, 1 17, 2941 

Tinney, C. G., Butler, R. R, Marcy, G. W., Jones, H. R. A., Penny, A. J., 

McCarthy, C, Carter, B. D., & Bond, J. 2003, ApJ, 587, 423 

Tinney, C. G., Buder, R. R, Marcy, G. W., Jones, H. R. A., Penny, A. J., Vogt, 

S. S., Apps, K., & Henry, G. W. 2001, ApJ, 551, 507 

Trilling, D. E., Lunine, J. 1. & Benz, W. 2002, A&A, 394, 241 

Udry, S., et al. 2000, A&A, 356, 590 

Udry, S., et al. 2003, A&A, 407, 679 

Vaughan, B. A., et al. 1994, ApJ, 435, 362 

Walker, G. A. H., Walker, A. R., Irwin, A. W., Larson, A. M., Yang, S. L. S., 
& Richardson, D. C. 1995, Icarus, 116, 359 
Ward, W.R. 1986, Icarus, 67, 164 
Ward, W.R. 1997, ApJ, 482, L211 
Yelle, R. V. 2004, Icarus, in press 



